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Abstract 



In this work, shells are mathematically constructed by applying the cut and paste 
procedure to D-dimensional spherically symmetric geometries. The weak energy con- 
_ dition for the matter on the shells is briefly analyzed. The dynamical evolution is 

^ ■ studied, and the general formalism for the stability of static solutions is presented. 

On . Several examples corresponding to diff'erent spacetime dimensions and values of the 

^ I parameters are considered. 

in 

O" 1 Introduction 
(N 

The study of thin shells in general relativity has been developed mainly in the framework of 
the formalism introduced by Darmois and Israel [1] . The central tool for analyzing the mat- 
^1 ter characterization and dynamics of surface layers are the Lanczos equations [HE], which 
■ relate the surface energy-momentum tensor of a shell with the jump of the extrinsic curva- 
ture tensor across it. Apart from cosmo logical applications, the Darmois-Israel formalism 
has been applied to highly symmetric configurations, as spherical and cylindrical shells. The 
linearized stability analysis of spherical shells was carried out by several authors (see [3lll] 
and the references included in these works). The formalism was applied to bubbles, shells 
around stars and black holes, and in the mathematical construction of traversable Lorentzian 
wormholes (see for example |SHZ] and references therein). Shells in more than four dimen- 
sional backgrounds have been considered in the construction of wormholes [2111] • Also, the 
Darmois-Israel formalism was applied to the collapse of a spherical dust shell into a Reissner- 
Nordstrom black hole in D spacetime dimensions |T0]. Here we address the characterization 
and general aspects of the dynamics of spherical shells in backgrounds of -D > 4 dimensions, 
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within the context of D- dimensional Einstein gravity and Maxwell electromagnetism. In 
section 2 we present the general formalism for spherical shells. In section 3 we briefly discuss 
the weak energy condition in relation with the masses and charges of the original metrics 
from which the construction starts. In section 4 we consider a non perturbative treatment 
of the dynamics and discuss some particular cases (an explicit solution is given in the Ap- 
pendix). In section 5 we present a general perturbative approach of the dynamics, suitable 
for a stability analysis of static solutions, and apply it to different spacetime dimensions; 
explicit results are shown and compared for four, five and six spacetime dimensions. Finally, 
the results are summarized in section 6. We adopt usual units such that c = 1. 



2 Mathematical construction 

We start the construction from two spherically symmetric D-dimensional manifolds Wl- and 
971+ with metrics 

dsl = -f±{r^)dtl + f^\T^)dTl + rld^l (1) 

(n = D—2) and boundaries S_ and E+. We identify the boundaries: S_ = S+ = S, and then 
paste the manifolds VJH- and TI+ at the hypersurface E, defined by J-'(r, r) = r — a(r) = 0, 
with T the proper time measured by an observer on the surface. The resulting manifold 
971 = 97t_ U 97t+ is geodesically complete and the corresponding line element is continuous 
across S as long as the coordinates in each side are related by f^{a)dt^_ = fj^{a)dt^. The 
induced metric on E is of course unique and has the form 

dsl = -dT^ + a'{T)dQl. (2) 

The joining of the two D-dimensional metrics implies a matter shell placed at r = a. As- 
sociated with this, we have a jump of the extrinsic curvature Kl across the surface, which 
is related to the energy-momentum tensor Si on the {D — l)-dimensional manifold by the 
Lanczos equations [2] 

(Kl) - K51 = -8nSi, (3) 

where (K^) = Kl_^_ — _, K = {6jKl). The components of the extrinsic curvature at both 
sides of the joining surface read 



K 



(4) 



- -^±m^ , (5) 

V f±{a) + a2 

where Ok {0 < k < n) are the angular coordinates, a prime stands for a derivative with 
respect to the radius and a dot means d/dr. The isotropy of the configuration leads to a 
simple diagonal surface energy-momentum tensor with energy density a = —S^ and pressure 
P = Sg^, given by: 

8na = ^V«' + /-H-^Va2 + /+(a), (6) 

fn-l\ ■d + f'_{a)/2 , a + /;(a)/2 
Sup = - Sncr , H , =, (7) 



2 



where we replaced a in the second equation for simphcity. These equations or any of them 
plus the conservation equation 

A(.O+P^^0 (8) 

are the starting point for the study of the shell dynamics (there is no problem with starting 
from a static metric of the embedding, because of the Birkhoff theorem). Note that as long 
as the condition /-(a) > /+(a) is fulfilled, the energy density a is positive. In what follows 
we will also explore in detail this point. 



3 Weak energy condition 

We will restrict our general analysis to shells of non exotic matter. Normal matter at 
the shell must fulfill the weak energy condition: a > and a + p > 0. If any of these 
inequalities is violated, the shell would be constituted by exotic matter. In this context, it 
can be convenient to define the auxiliary quantity f2„ = 27r^"+-^''/^/r[(n+ l)/2]. Thus, within 
the Einstein-Maxwell framework, the general form of the function /(r) for a black hole in 
D = n + 2 dimensions is (see [11] ) 

f(r) = l-^ + ^-^ (9) 

■1^' ) ^ ^n-l ^ ^2(n-l) 3 V^'' 

where A4 = 8TTGDiTi/{nQn) and = 2q'^/[n{n — 1)] with m and q the ADM mass and 
charge, and A is the cosmological constant; Go is the Newton constant in D dimensions. 
Therefore, for static solutions (a = constant) the weak energy condition to be satisfied by 
the energy density and the pressure leads to 



STia = '^^j4a)--^U{a)>Q, (10) 

(X 



Stt cr + p = — 7ra H , > 0- 11 

Note that while in the case q± = ^ and A = the weak energy condition is trivially 
satisfied if m+ > m_, this is not the situation in general. Normal matter requires at least 
/-(a) > /+(a). Examples in which the weak energy condition would be violated are easy to 
find: Consider a charged bubble (m_ = 0, g_ = 0, m_|_ > 0, g+ 7^ 0) in a five dimensional 
background with vanishing cosmological constant. The condition to be outside the horizon 
is 8G'5m+/(37r) < + (3'^/(3a^); but this range includes 8G5m+/(37r) < g^/(3a^) (note that 
this does not imply a naked singularity, because here we are not speaking about a point 
mass, but about a shell of finite radius). Now, if this is the case we have /+(a) > 1 = /_(a), 
which means that the energy density would be negative if < nq'^/lSG^rri-^-). Then for 
given charges and masses there is a lower bound for the possible bubble radius. 

A similar care should be taken in more general cases: if the cosmological constant is 
(reasonably) assumed to be equal inside and outside the shell, one should always start the 
analysis from a static solution satisfying the necessary condition 

2M- - QW''' < 2M+ - QW^"". (12) 
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Besides, one should also ensure (T + p> 0. Now, from the conservation equation ([8]) it is not 
difficult to show that 

a + p = —aa' /n; (13) 

in practice, then, for a non trivial form of f±, the most simple way to deal with the issue of 
avoiding exotic matter is to numerically find the intersection of both conditions a > and 
a' < 0. This will be the procedure followed in the stability analysis of section 5. 



4 Nonperturbative approach: examples 



The dynamics of the shell is straightforwardly obtained from the last equations of section 2. 
Squaring twice and rearranging expression Q, we have 



+ V{a) = 0, 



(14) 



which has the form of the energy conservation of a point particle in a one dimensional 
problem, with the function V{a) playing the role of a "potential". This potential has the 
form 



V{a) 



/_(«) + /+ (a) 



n(/_(a)-/+(a)) 



- 2 


AiTaa{a) 




n 



(15) 



16'Kaa{a) 

In order to obtain an equation of motion for a, the dependence of a with the shell radius 
should be given. The conservation equation ([8j) leads to aa' +n{a+p) = 0. Once an equation 
of state yielding a relation p = p{a) is adopted, we can integrate to obtain 



In a 



da 



a 



(16) 



This integral should be inverted to have a (a). From Eq. (IT^ one has da/dr = =Pa/— l^(a) 
(recall that from Eq. f[T^ the function V{a) must be negative along the allowed range of 
radii), which once a{a) is known can be integrated giving 



da 



+ C2. 



(17) 



Then this relation must be inverted to obtain the solution a(r). In brief: for a given equation 
of state, Eqs. (IT^ . ( IT^ and ( 1T7|) give the general solution of the problem. Clearly, for a 
general metric function most equations of state will lead to an analytically non tractable 
problem. However a qualitative analysis is possible for some physically meaningful exam- 
ples; in the rest of this section we will restrict the discussion to the case of five spacetime 
dimensions. 

1. Case A = 0, non charged shell of dust or non relativistic matter {p <^ a) around a 



charged object (g+ 
the potential reads 



g_ = q). In this situation a 



We set M = 47raV/3 and 



V{a) 



47^2 



^(j5(m_ — 
3M 



a4 V 3 
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2. 



Assuming that both m_ and m+ are positive, then we have a monotonous potential 
for g^/S < Af^, and a potential with a minimum for > M^. In the first case an 
initially static shell can only collapse, while in the second case an oscillatory motion 
could, in principle, take place. Note that if g = no oscillations are possible. 

Case A = 0, non charged shell around a non charged black hole (g_ = g+ = 0), linear 
equation of state: p — rja;we assume < 77 < 1 in order to allow for the interpretation 
of 77 as the squared velocity of sound on the shell. In this situation the conservation 
equation leads to the behavior of the energy density a{a) ~ ka~^^^~^'^\ and the potential 
has the form 



V{a) 



8G'5(///_ + ///+) 1 



An 



(19) 



Because we assume < 77 < 1, now we have a positive power of the radius; however, 
all the powers appear with a negative coefficient, so that an unbounded motion is not 
forbidden. 

3. A charged bubble (m_ = g_ = 0, m+ = m, q'+ = q) of dust or nonrelativistic matter 
{p <^ a). With the definition for M introduced above, the potential reads 



V{a) = 1 



SGsm 
37ra2 



6a4 



47r2 



SGgm - 7rgV(2a2) 



3M 



M2 

-F- (20) 



If the bubble has vanishing charge, the collapse is unavoidable, because the potential 
is a monotonically increasing function of the shell radius, and V{a) goes to —00 when 
a goes to zero. For a non vanishing charge, the presence of the associated positive term 
in the potential could allow a different behavior. 

4. A dust shell {p = 0) around a black hole, in a cosmological constant background 
(A 7^ 0). The potential for such a model has the form 



V{a) 



1 - 



8G5{m-+m+) Aa^ 



47r2 



8^5(771- — 77T.+ ) 
3M 



M2 



(21) 



Here the interesting situation is that of a negative cosmological constant, because 
for A < (anti-De Sitter background, no cosmological horizon present) the potential 
diverges when a — > 00, and an unbounded evolution is then excluded. 

In the Appendix we will consider another example for which we will give an analytical so- 
lution for the dynamics under certain approximations. In the following section, instead, we 
will present a general procedure for studying the stability of static solutions under pertur- 
bations preserving the spherical symmetry. The approach will have the nice feature of being 
independent of the equation of state for the matter on the shell. 
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5 Perturbative treatment: stability 



As a less general but physically sound approach to the dynamics of the shell we can consider 
small perturbations preserving the symmetry around a static solution. Our procedure will 
then be similar to the treatment in Refs. [3111] • As noted in the preceding section, once an 
equation of state is adopted we can formally take the potential as a function of the shell 
radius. For a perturbative treatment of the stability of static solutions it is enough with 
the analysis of the first and second derivatives of V{a) at a radius for which d = 0. 
Equilibrium satisfies V{aQ) = and V^'(ao) = 0, and stability requires V"{ao) > 0. We define 
the functions 

S{a) = f-^-^ + f^^-\ R^a) = f-^-^'^f^^-^ (22) 

and we introduce the relation M = a"f2„(T, with M not necessarily a constant but a function 
of the radius. By defining u = 'n,fi„/(47r) the subsequent expressions can take a simpler 
form. Then the potential can be written as 

The first and second derivatives read 



Via) 



S" 



2 f M 



u;2 V a"-i 



/2 



2M / M 



2r,n-l \ r,n-l 



71—1 / „n—l 



-2u'^R'R 



M V M 



a" 

uj" [ a 
~2 



co^R^ 



M 



l2 



M \ M 



M 



R'^ + RR" 



(25) 



In the equations that follow it must be understood that the functions are evaluated at an 
equilibrium radius oq and the prime means the derivative with respect to ao. Equilibrium 
implies V(ao) = 0, from what we can express 



M \' 



2 „n-l 



u a, 







2M 



, ui^R'R f Aq 



n-l\2 ,,2p2„n-l 

^ u rC a, ' " 



M 



2M 



M 



X(ao). (26) 



Stability implies V"{a^ > 0, which with this definition leads to 



2M f M 



where 
y(ao) = S" 







+ 



2M 



M 



M 



n— 1 / „n— 1 



<r(ao)-^X2(ao), 



2 D' D "0 



2uj^R'R 



M V M 



bj I a, 
~2 



(27) 







M 



R'^ + RR" 



(28) 
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Second derivatives of M would lead to second derivatives of the energy density. But with 
the aid of the conservation equation we can write 



ao 



[n + l)cr' + np'] . 



(29) 



This seems to introduce the pressure, and consequently the choice on an equation of state, 
in the formalism. However, because the expressions above apply to the equilibrium, we can 
define the parameter r] = p'/cr', with both derivatives evaluated at the cquihbrium radius, 
so that only the lowest order of the expansion of the equation of state appears. Then going 
back from cr to M we obtain 



M 

^n— 1 



n 



ao 



M 

^11—1 



M 



^0 J 



1 + 



n 



n 



M 



M 



M 



l2 



+ ■ 



n 



,n— 1 



M 



M 



M' 



Introducing the definition 



Z{ao) = u^R^ 



,n— 1 



M 



n-1 

ao 



M 

1 



n-1 



M 



2M 







M 

,n— 1 



M 



the condition for stable equilibrium can finally be put in the form 



1 + 



n 



n 



(30) 

(31) 



1 + 



n 



n 



1 



77 (32) 



Z{ao) < Y{ao) - ^X^{ao). 



(33) 



The subsequent analysis needs in general be carried out numerically, and can be performed 
in terms of the parameter 77, which can be understood as the square of the velocity of sound 
on the shell. Then, the preferred range would be < 77 < 1. 

As an application of the formalism, we have studied the stability of charged bubbles 
(r?T._ = 0, 5_ = 0) and charged shells around non charged black holes (r?T._ = 0.5r?7,_|_, — 0) 
in four, five and six spacetime dimensions. We have set the shell radius beyond the horizon 
radius of the original outer manifold (so this horizon is removed); which is given by 



l/{n-l) 



(34) 



where Ci = STi/inVLn) and C2 = 2/[n{n — 1)]. We have restricted the analysis to shells of 
normal matter, i.e. matter satisfying the weak energy condition. The results are displayed 
in Figs. 1-6. We have not restricted the graphics to < 77 < 1, though the results within 
this range are of more physical interest. We can see that in all cases the stability regions 
change their shape when the charge reaches the critical value qc = Gom+Ci/ \fc2. from which 
the horizon in the original outer manifold disappears. For charges below qc-, the regions of 
stabihty for bubbles appear to be shghtly larger than those for shells around black holes. 
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l9+l=0 



19+1=0.8 q. 



l9+l=0 



19+1=0.8 q. 



\q+\=qc 




\q+\ = \A q. 




12 3 

<!o(G4'"+)"' 



Figure 1: Stability regions (grey) for 4-D 
charged shells around vacuum; Qc = G^rrij^. 



19+1=0 



19+1=0.8 q. 





l<?+l = l.l 9t 




Figure 3: Stability regions (grey) for 
5-D charged shells around vacuum; 
qc = 4G5m+/(7rV3). 



2 3 4 



2 3 4 

00(04014)"' 

l9+l=l.l 9c 




ao(G4m+) ' 



Figure 2: Stability regions (grey) for 4-D 
charged shells around a non charged black 
hole; Qc = G^m^. 



l9+l=0 



19+1=0.8 q. 




Figure 4: Stability regions (grey) for 5-D 
charged shells around a non charged black 
hole; Qc = 4G5m+/(7r\/3). 



while the reverse takes place for charges beyond Qc, however, for q > Qc the stability regions 
begin at a slightly larger radius as the charge increases. For any number of dimensions, 
stable bubbles with r) < require q > q^, while stable layers with 77 < around black holes 
are possible for q > qc- The most interesting feature of the results is that the stability 
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l9+l=0 



19+1=0.8 qc 





ao{G(,m+) 




lg+l=o 



l<?+l=0.8 q. 



Figure 5: Stability regions (grey) for 
6-D charged shells around vacuum; 
Qc = 3^/3Gem+/{2^T^/2). 
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\q+\=qc 



ao(Gsm*r"-' 



k/+l=i-i * 



Figure 6: Stability regions (grey) for 6-D 
charged shells around a non charged black 
hole; = 3V3GQm+/{2nV2). 



regions for fixed |g+|/gc are larger as the number of dimensions increases. In particular, 
as the dimensionality of spacetime increases, for fixed ao{G D'm+)~^^^"'~^^ stability becomes 
compatible with smaller values of the parameter r]. 

The procedure above is also valid for wormholes if the outer part of both manifolds is 
taken, and the signs in the expressions of the energy density (E]) and pressure ([7]) are suitably 
changed. In particular, if the two manifolds from which the construction starts are equal 
copies of a five dimensional spherically symmetric geometry, results including those of Ref. [8] 
can be recovered. In this sense, this section constitutes a generalization of the procedure 
presented before in that work. 



6 Summary 

We have addressed general aspects of the characterization and dynamics of spherically sym- 
metric shells within the framework of Einstein gravity and Maxwell electrodynamics in 
D-dimensional spacetime. We have applied the Darmois-Israel formalism extended to D 
dimensions to the mathematical construction of the shells, starting from spherically sym- 
metric geometries associated to, in general, charged black holes with a cosmological constant 
background. We have discussed the conditions to be imposed on the parameters in order 
to ensure that the matter on static shells is not exotic, i.e. it satisfies the weak energy 
condition. Then we have considered the full dynamics preserving the spherical symmetry; 
we have given the formal general solution of the problem, and discussed some examples by 
means of the analogy of the equations of motion and the potential of a point particle. Also, 
an explicit solution is given in the Appendix. Finally, we have presented a formal approach 
to the study of the stabihty of static shells under spherically symmetric perturbations; this 
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is given in a concise fashion suitable for the actual application, which in general needs to be 
carried out in a numerical form. We have applied the formalism to obtain explicit results 
regarding the stability of charged bubbles and shells around black holes in four, five and 
six spacetime dimensions. The most interesting result seems to be that as the number of 
spacetime dimensions increases the stability regions in parameter space become larger. 



Appendix 

Consider the five dimensional metrics with f±{r) = l — 2m±/ (vrr^), in which we have adopted 
units such as G5 = 3/4. We assume p = 0, corresponding to the physically interesting case 
of a shell of nonrelativistic matter or dust. This assumption implies an energy density with 
the dependence 

k 



This leads to the equation of motion 



.2 m_. + m+ f — m+\ (Ank/S) 



(35) 



The time evolution of the shell can be obtained by integrating Eq. (136 p . which gives the 
proper time in terms of the shell radius. Defining a = — 1 + 9(m_ — m+)^/(87r^A;)^ and 
f3 = —{rri- + m+)/7r we have 

r = ±! , +C. (37) 

J ^aa'^ - Pa^ + {A-Kk/?,Y 

Inverting the relation given by the solution of this integral yields the shell radius as a function 
of the proper time. The integral can be expressed in terms of elliptic functions; however, a 
qualitative analysis can be carried out by recalling that the problem is formally analogous 
to the dynamics of a point particle in a potential 

(47rA;/3)2 , , 

with a null total energy. For /3 < only an accelerated contraction or a decelerated expansion 
are possible; however, for /3 > the possibility of an unbounded accelerated motion exists, 
for suitable values of the parameters. It is easy to show that for /3 > the potential has a 
maximum for 

4 [2 

Om = -'^kJ -, (39) 

so that V'{a) < for a > a^. Values of the parameters m_ and m+ can be chosen so 
that a > 0, /9 > 0. In this case, depending on the relation of a and /3 with k > 0, V{am) 
can be both positive or negative, while the potential asymptotically tends to the constant 
—a < (a simple numerical analysis shows that, taking a scale such that m+ = 1, for 
example in the case m_ = —2 and k = 0.11 we have a positive maximum, while for k = 0.1 
the maximum is negative; it is also easy to see that with such values of the parameters 
the condition am > A/2m+/7r is fulfilled). Thus, for such sets of parameters and an initial 
condition d, > 0, we have three kinds of evolution: 
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1. If y(am) is negative, for any initial condition Oj < (always beyond the radius 
■\/2m+/7r) the shell radius undergoes a decelerated increase until a = Om is reached, 
and beyond this point the subsequent evolution is accelerated. 



2. If l^(ain) is negative but > Om, the shell expands with positive acceleration. 

3. If, instead, V{a-[a) is positive, the potential has two positive real roots, and for above 
the largest root the shell undergoes an accelerated expansion. 

Because the derivative of the "potential" vanishes for a — > oo, in all cases the positive 
acceleration decreases as the radius grows. In brief, the model undergoes an accelerated 
expansion with decreasing acceleration. This kind of evolution of the shell radius is possible 
only if /3 > 0, which implies a negative mass. But this condition can be fulfilled without 
any exotic matter in the (3 + l)-dimensional manifold resulting from the cut and paste 
construction: if m_ < 0, m+ > and |m_| > then cr > and also cr + p > 0. The 

exotic matter -negative mass m_- is placed under the radius a in the (4 + 1) -dimensional 
manifold from which we started the construction. Therefore a (3 + l)-dimensional shell 
with matter in the form of dust of positive energy density is compatible with an accelerated 
expansion. Recall that given the form of the induced (3 + l)-dimensional metric one could 
understand the result for the shell as the time evolution of a closed cosmology which, quite 
interestingly at present, undergoes an accelerated expansion. This accelerated toy "universe" 
would present two positive features: the unusual kind of matter driving the acceleration is 
not placed within (3 + l)-dimensional spacetime, and the framework in which such kind of 
evolution is possible is that of classical general relativity. Of course, this is not to be taken 
too seriously; i.e. we do not pretend that this observation has a phenomenological interest. 

An approximate explicit solution for the time evolution of the model can be obtained for 
the regime in which the radius a is large enough. More precisely, for o; > and as long as 
aa^ — l3o? :$> (47rA;/3)^, we can expand 



-1/2 



2(2^/,-/3)^ 



(40) 



Under this approximation the integral form of the time as a function of the radius is 



a da / 2 , 

^/aa? - (3 V3 



The integrals are easily calculated, and the result is 



a^da 



(aa4 - /3a2)3/2 



+ C. 



(41) 



a 



3/2 



P^Ja? - P/a \/3j 



+ I — arccos 



+ C. (42) 



The approximate evolution of the shell radius is obtained inverting this relation. Note that 
for very large values of a the evolution is almost linear with time (the zeroth order result 
is given by the first term, and is simply r ~ y^a"^ ~ P/a), which is consistent with the fact 
that, as pointed out above, the derivative of the potential tends to zero when a tends to 
infinity. 
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